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Abstract. We study the field theoretical model of a scalar field in presence of spacial 
inhomogeneities in the form of one and two finite width mirrors (material slabs). The 
interaction of the scalar field with the defect is described with position-dependent mass 
term. 

For a single layer system we develop a rigorous calculation method and derive 
explicitly the propagator of the theory, the S-matrix elements and the Casimir self- 
energy of the slab. Detailed investigation of particular limits of self-energy is presented, 
and connection to known cases is discussed. 

The calculation method is found applicable to the two mirrors case as well. With 
its help we derive the corresponding Casimir energy and analyze it. For particular 
values of parameters of the model an obtained result recovers the Lifshitz formula. 
We also propose a procedure to unambiguously obtain the finite Casimir self -energy 
of a single slab without reference to any renormalization conditions. Wc hope that 
our approach can be applied to calculation of Casimir self-energies in other demanded 
cases (such as dielectric ball, etc.) 



PACS numbers: 03.70.+k, ll.10.-z 



1. Introduction 



Usual models of interaction of elementary particles are constructed in a framework of 
the quantum field theory (QFT) in the homogenous infinite space-time pp. However, as 
such they can not be applied for description of phenomena of interaction of quantum 
fields with macroscopic bodies. At least the form of latter must be presented in the 
model, and it can change essentially both the spectrum and dynamics of the excited 
states of the model as well as the properties of the ground (vacuum) state of the system. 

First quantitative results for such effects were obtained by H. Casimir in 1948. He 
predicted [2] macroscopical attractive force between two uncharged conducting plates 
placed in vacuum. The force appears due to the influence of the boundary conditions 
on the electromagnetic quantum vacuum fluctuations. Nowadays the Casimir effect is 
verified by experiments with a precision of 0.5% (see [7j for a review). 
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The properties of vacuum fluctuations in curved spaces, the scalar field models 
with various boundary conditions and their application to the description of real 
electromagnetic effects were actively studied throughout the last decades (see discussion 
and references in [I]- [9]). However, it was well understood from the beginning that 
boundary conditions must be considered just as an approximate description of a complex 
interaction of quantum fields with the matter. A generalization of the boundary 
conditions method has been proposed by Symanzik [3]. In the framework of path 
integral formalism he showed that the presence of material boundaries (two dimensional 
defects) in the system can be modeled with a surface term added to the action functional. 
Such singular potentials with 5-type profile functions concentrated on the defect surface 
reproduce some simple boundary conditions (namely Dirichlet and Neumann ones) 
in the strong coupling limit. The additional action of the defect should not violate 
basic principles of the bulk model such as gauge invariance (if applicable), locality and 
renormalizability. 

The QFT systems with 5-type potentials are mostly investigated for scalar fields, 
see for instance [10]. In [II] [13] the Symanzik's approach was used to describe similar 
problems in complete quantum electrodynamics (QED), and all 5-potentials consistent 
with the QED basic principles were constructed. 

It seems quite natural to try applying the same method for the description of the 
interaction of quantum fields with bulk macroscopic inhomogeneities (slabs, finite width 
mirrors, etc.) and to study Casimir effects in systems of such kind. The Symanzik's 
method of adding extra potential terms into the action of a system was used to model 
the interaction of quantum scalar fields with bulk defects in a number of papers (e.g. 
[26] . [29]- [32], and others). However, most of them were devoted to study of a limiting 
procedure of transition from a bulk potential of the defect to the surface (^-potential 
one, without paying much attention to other properties. 

Traditionally, there are several approaches to study electromagnetic Casimir effect 
between material slabs. The electromagnetic field and the material bodies can be 
treated macroscopically and by employing a dissipation fluctuation theorem one obtains 
the field correlation functions that are needed to construct the Maxwell stress tensor. 
In alternative approach the presence of dielectric bodies is described by means of a 
spatially varying permittivity that is a complex function of the frequency. Within these 
approaches the influence of dielectric and geometrical properties on the Casimir force 
between dispersing and absorbing plates has been studied in j22j El]. A dependence of 
the Casimir force on thickness of the slabs was considered in [25]. A calculation of the 
Casimir force in a (multi)layer system have been also performed in [T9 | |20~ 1 |2T | [23] . Most 
recently there was an attempt to construct quantum electrodynamics in the presence 
of dielectric media (i.e. volume defects of special kind) [18]. From the field theoretical 
point of view none of these methods were truly successful. Particularly the celebrated 
Lifshitz approach [TJ] to the description of the interaction between dielectric bodies still 
attracts a lot of arguments if dispersion is present in the system and formally is proved 
only for the interaction over the vacuum gap, see [15]. While in some approaches there 
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could be no action functional constructed [5J, there are also general severe problems 
in construction of path-integral quantization for the QED systems with dielectrics [16] 
and calculation of the heat kernel coefficients for such systems [28J. Moreover, as yet 
remains unresolved the mystery of Casimir energy of a dielectric ball which is defined 
unambiguously only in the dilute limit or for the speed of light being constant across 
the boundary [T7] . 

On the other hand, existing results for the Casimir energy of a scalar field in 
presence of a single planar layer of finite width £ are contradictory. The formulae 
presented recently in [33] do not coincide with previous calculations made in [26] as 
discussed in [33]. The only attempt to calculate the propagator in such system was 
undertaken in [35], where hardly any explicit formulae were after all presented. The 
system of two interacting slabs has never been considered within QFT models. 

Thereby, one can see that the specificity of finite volume effects generated by 
inhomogeneities in QFT has not been yet adequately explored. Our work is dedicated 
to clarify the problem, and to solve existing controversy within an accurate and 
unambiguous approach. On the other hand, considered type of square well potentials 
has its own range of applicability including transport in graphene and other condensed 
matter systems [39], and (apparently) superluminal effects in wave propagation [40J. 
While consideration of the scalar fields may only appear restricting, for the case of 
planar geometries it is apparently sufficient. It is due to the well known fact that in this 
case the theory of electromagnetic (EM) field is equivalent to two scalar field theories 
as TE and TM modes of the former field do not mix. The study of EM fields in the 
geometries, when such decomposition does not appear, will be the scope of the future 
work. 

Thus, we consider a model of massive scalar field interacting with the volume defects 
- finite width material slabs or mirrors, — whose properties are effectively described 
by single macroscopical parameter, the coupling constant. By construction the model 
is renormalizable in the sense applicable to QFT with spatial inhomogeneities given 
in [3]. Working with renormalizable models has an important advantage that their 
long distance (macroscopical) properties such as Casimir interaction can be described 
without explicit knowledge of microscopical structure even if there is a true physical 
cut-off present in the system, [4T] . Such description can be given in terms of just a few 
effective parameters. In our case, the role of these parameters is played by the coupling 
constant. Any physical cut-off dependence is disregarded in light of renormalizability of 
the model. 

We develop a purely field theoretical approach and construct explicitly the 
propagator for the system with a single finite width mirror (Sect. 13.11) and the S-matrix 
elements (Sect. 13. 2p . Based on these new results in Section |3~U| we recover the known 
expression for the self-interaction of a single slab and supply it with the investigation of 
its asymptotical behavior in different physically interesting regimes (Sections I3.5H3.7]) . 

In Section H] the Casimir energy of two interacting finite width slabs described 
within the scalar field model with step potential is obtained for the first time (Sect. 14.21) 
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and several limiting cases are considered (Sect. 14 . 3 p . We also reveal its correspondence 
to other known results in the literature, and in section 14.41 discuss its connection to 
the Lifshitz formula. Possible ways to construct divergence free limits for Casimir self- 
interaction of a solitary body are presented in Sect. 14.51 In the Conclusion the main 



results are summarized. In the Appendix A we briefly discuss a by-side question of 



comparison of different regularizations within QFT. 
2. Statement of the problem 

Let us consider a model of a scalar field interacting with a space defect with nonzero 
volume. Using the Symanzik's approach, we describe such system by the action 
functional with an additional mass term being non-zero only inside the defect 

S = S + Sdef (1) 

So = lJ d A x(j)(x)(-dl + m 2 )(f)(x) 

Sdef = \ j d 4 x6(£,x 3 )(j) 2 (x) 

where d 2 = d 2 /dxl + . . . + d 2 / dx 2 ^) . In the simplest case, the defect could be considered 
as a homogenous and isotropic infinite plane layer of the thickness £ placed in the x\x% 
plane. In this case it is sometimes called 'piecewise constant potential'. The distribution 
function 8(£, Xs) is then equal to l/£ when \x 3 \ < i/2, and is zero otherwise. In terms 
of the Heaviside step-function it can be written as 

9(1, x 3 ) = [6(x 3 + £/2) - 9(x 3 - £/2)}/£. (2) 

In the framework of QFT such potentials were approached for the first time in [26], and 
later in [30]- [33]. However, the latest consideration |33j was proved to be inconsistent 

To describe all physical properties of the systems it is sufficient to calculate the 
generating functional for the Green's functions 

G[J\ = N J D<p exp{-S[0] + J0}, A^ 1 = J D(f) exp{-S [^]} (3) 

where J is an external source, and normalization for the generating functional has been 
chosen in such a way that G [0] | a=o = 1 • 

Introducing in ([3]) an auxiliary field ip defined in the volume of the defect only, we 
can present the defect contribution to G[J] as 

exp { I dx J ^ dx^4> 2 (x) 1 = (4) 



C J Dtp exp I J dx J dx 3 ^— %^ + iyfnijjipj 



% We operate in Euclidian version of the theory which appears to be more convenient for the 
calculations. 
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which is written for the case of a single layer. C is an appropriate normalization constant, 
and k = X/£, here and below we denote with arrows the (0, 1, 2) components of 4- vectors, 
i.e. x = (xo,xi,X2). When generalized to other defects the ^-integration in (jlj should 
be performed only within their support. 

With help of projector onto the volume of the layer O = 9(xs + £/2) — 6{x^ — £/2) 

£ /2 

acting as ipO(J) = J dx f_ e , 2 dx 3 tjj(f), we can perform the functional integration over 0, 
and consequently over ip. As the result we get 

G[J] = [DetQ]- 1/2 e5 J ^ J , S L = D - k{DO)Q~\OD), (5) 
Q = 1 + k{ODO). (6) 

Here the unity operator 1, as well as the whole Q, is defined in the volume of the defect 
only [—1/2,1/2) x M 3 , and D = {—d 2 + m 2 ) -1 is the standard (Feynman) propagator of 
the free scalar field. We shall note here that the outlook of ()5]) completely coincides with 
the expression for the generating functional G[J] in the case of delta-potential defect 
instead of patchwise constant one subject to appropriate redefinition of the projecting 
operator O. It is also evident that a straightforward generalization is possible for non- 
constant k(— X/£) with A depending on x$. 



3. One mirror system 

3.1. Calculation of the propagator 

To calculate the propagator S defined according to (JSJ) let us first derive an explicit 
formula for the operator W = Q~ l . 

For this purpose we make the Fourier transformation on the coordinates parallel to 
the defect (i.e. xq, x\, x 2 ). In such mixed p-x 3 representation the propagator D of the 
system without a defect is given by 

f d 3 n ~ P -y/V\x\ 
D(x) = J -^e^V^M, V v (x) = — (7) 

with E = \/p 2 + m 2 , p = (po,Pi,P2)- Then W can be defined through the following 
operator equation 

W + kV E 2W = 1. (8) 

By construction the free scalar propagator V E 2 (x, y) = V E 2 (x—y) is the Green's function 
of the following ordinary differential operator 

K v (x,y)=(~ + v\6(x-y) (9) 

for V = E 2 . Multiplying both sides of (jSJ) with K E 2 and using obvious relation 
K v - Kyi = V - V we get 

K P U = — K (10) 

where p = k + E 2 and U = W — 1 . 
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The general solution to this (inhomogeneous) operator equation can be written as a 
sum of its partial solution and the general solution of its homogeneous version. Taking 
into account the symmetry property of U: U(x,y) = U(y,x), which follows from its 
definition, we arrive at 



U(x, y) = -kV p {x, y)+ae {x+y) ^~ p + b (e^)^ + 6 (v-*)Vp) + ce -0 B +i'Vp(ii 

where a, b and c are some constants. To derive them we introduce W = 1 + U into (E 
and require its identical validity for all x and y. It yields 

^ 2 e^ , - Vk + E 2 ) 2 



2\/k + E 2 ' 2-s/k + E 2 

1 



(12) 



S(P, x 3 ,y 3 ) 



e 2i^+Wf E + ^ k + e 2 ) 2 — (E — Vk + E 2 ) 2 ' 

Substituting Q -1 = 1 + U into (jSJ) and using (fl2l) we can finally derive the explicit 
formulae for the modified propagator of the system S = S(p,X2, y 3 ) 

' S^, x 3 < -£/2, y 3 < -£/2 

S_ D , x 3 <-£/2, y 3 e{-£/2,£/2) 

S_ + , x 3 < -1/2, y 3 >£/2 { > 

Soo, x 3 e (-£/2,£/2), y 3 e (~£/2,£/2) 

where 

= e~ E ^-y^ E{e+X3+m) _ (14) 

2E 2E v ; v ' 

S _ o = £ (V^+^^+^v^T^ -E) + e { - u ' 2 -^ V ^{^+E 2 ' + £)) e E ^ + ^ (15) 

S_+ = 2£V^TE 2 e ( ^+ W+E ^ +E( > X3 -^ (16) 



Soo = . j 2/ t cosh[(x 3 + i/3)v / ^ + ^]+ (17) 

2v K + £/ 



+ ( v ^T^-^) 2 e^^(l"l-^ + (E + v / ^T^) 2 e v/ ^ ! '- --^ ! 

We have divided the general expression of the propagator into four parts according to 
the position of x 3 , y 3 relative to the defect, and all other cases could be easily derived 
using the symmetry properties of the propagator. 

Despite that the problem of square-well potential has been studied extensively in 
quantum mechanics, its field theoretical interpretation was lacking. To the best of 
our knowledge the only attempt to calculate the full propagator for such system was 
presented in [35]. However, the final explicit expression is lacking there, but instructions 
for its construction are given to the reader. Following these instructions one can compare 
[35] with (Tl4l) -(!T7l) and find them coinciding up to the Wick rotation. One should also 
note that in Euclidian version of the problem the difference between under- and above- 
barrier scattering vanishes. The 5 00 , for 2-dimensional problem were also presented 

in [8]. 
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Fixing position of one of the interfaces of the potential and taking the infinite 
width limit transforms the potential into a step-function one. The propagator for non- 
relativistic problem with such potential was calculated earlier, e.g. [37]. Subject to 
necessary redefinitions of parameters, our results (1T4T) - (TTTT) for £ — > oo coincide with 
ones presented there. 



3.2. Scattering on the slab 

In the quantum field theory the scattering processes are described by the S-matrix 
elements between different states. To construct them for one-particle states of particular 
impulses, we consider the full propagator ([13]) in real (Minkowski) space, and cut the 
external lines of the free field propagator 

H(x, y) = D-\x, z')S{z', z")D-\z" , y). (18) 

The defect potential depends only on x 3 thus both the energy (p ) and impulse parallel 
to the defect (p) are conserved. We omit corresponding delta-functions in the Fourier 
representation. In the transverse direction we have 



H(p 3 ,q 3 )= J dx 3 dy 3 e tp ^ +l ^H(x 3 ,y 3 ) 



= - k f' 2 dx 3 f' 2 r///, (6(x 3 - y 3 ) + U M (x 3 ,y 3 )) 

J-e/2 J-e/2 

where the last equation is due to the definition of S and the obvious property that 
Ugl(x, z)V E 2(z, y)0(y) = S(x — y)0(y) is only non-zero when both its arguments lie 
in the defect region (x,y G (-1/2, £/2)). By U M we denote the U operator ([11]), (fl2j) 
subject to the inverse Wick rotation. 

Further we note, that in the Minkowski space the on-shell condition leaves only two 
possible values for transversal impulse p 3 : p 3 = ±E, E = \/pq — p 2 — ni? . Due to the 
symmetry reasons there are only two independent values of H depending on the sign 
°f P 3 , Q 3 — iE . Thus the non-trivial S- matrix elements depend only on one parameter 
p = E and are given by 



8p 2 Jn - p2 e ^p+V^) 

hj — — *2zp — (19) 

( e 2V^7 - l)( K - 2p 2 ) + 2ipv^T^2(l + e^v 7 ^) 

fc = — ^(l-e^)e^ 

(e^V^ - 1)( K - 2p 2 ) + 2i Py /K-p 2 (l + e 2e V*~^) 

where hf = H(p,—p) is the forward scattering amplitude, and h b = H(p,p) — the 
reflection one. 

It is possible to check that the presented S'-matrix amplitudes satisfy the following 
analog of the optical theorem 

2 lm{h f ) = -^(IM 2 + N 2 ), 2 MM = -^(hfK + h}h b ) (21) 
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3. 3. Casimir Energy 

The Casimir energy density per unit area of the defect S can be presented with the 
following relation 

£ = -^lnG[0] = ^Trln[Q(x,y)]. (22) 

In the second equality we used (EJ), T is the (infinite) time interval and S — the surface 
area of the defect. For the explicit calculations we first make the Fourier transformation 
as in (J7|). Then 

£ = A d J ^jliMQ^^lft,)], (23) 

where we introduced dimensional regularization and an auxiliary mass parameter fi to 
handle the UV-divergencies. 

Using the definitions of U = U(k) and Q = Q(k) we can express the K-derivative 
of the integrand of (123!) as d K \n.Q = T> E iW = —Uk" 1 . Then for the energy density we 
get 

£ = -^f^/50^' (24) 

We have chosen the lower limit of integration over k' to satisfy the energy normalization 
condition £\ K =$ = 0. As we show below the integral is convergent at k' = 0. 
The trace of the integral operator U(k) is straightforward 

f^ 2 , s 4a sinh(£,/p) -Ik 
TrU= dx U(x, x) = 2bi + )-^- (25) 

J -1/2 



where we already used that a = c. Using a and b given in (1121 . one easily notes that 
TrU ~ —£k/(2E) when k — > 0, thus supporting the above statement. 
Next, putting fjX2 jl into ([25]) we can prove directly that 

- e -i(E+^p)- 



TrU = —K-— In 



(26) 



dK [ AE^pi 

Thus, from (|24|) and ( 126]) we obtain the following expression for the Casimir energy 

s = "" d I £^ in K (eCV ~" (E + " e ^ iE " ^ )2) 1 (27) 

For the first time the explicit results for the Casimir energy of single slab were obtained 
in [26], and recently rederived in [27J. Upon using different regularization scheme and 
notation they coincide with <\27L The question of comparison of results obtained in 



different regularization schemes is addressed in Appendix A 



To extract the UV divergencies at d — 4, we break up the energy in two parts 

£ — £fin + £ div 5 (28) 

where 

1 f°° 

£fin = 4^ E ^P 2d P> ( 29 ) 
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H(p) = In 



A 



4££ 2 



A/z 



4-d 



2(2 



7T 



id— 1 



d d ~ l p 
2E 



1 



A 



A£E 2 



(30) 



It easy to check that the first item, Efi n , is finite if we remove regularization, while Ediv 
is divergent. Its asymptotic in the limit of removal of regularization is given by 

A(2£m 2 + A) 



E 



dvu 



32n 2 £(d - 4) 
A 



+ 



Q4tt 2 £ 



A + (2£m 2 + A) 7b - 1 - ln(4vr) + 2 In 



m 



(31) 
+ 0(d-4) 



Here 7^ = 0, 577215 is the Euler constant. Explicit calculation of Ediv will be used in 
the considerations of Section 13.71 



S.Jf.. Renormalization procedure 

Simple analysis of the dependence of Edw ( l3"Uj) on the parameters A, £ shows that for the 
renormalization of the model at the one-loop level considered here we must add to the 
action at least the following field-independent counter-term SS 

5S = f\ + g\ 2 r\ (32) 

with bare parameters / and g (of the mass dimensions two and zero correspondingly) 
[§|. It allows us to choose these parameters in such a way that the renormalized Casimir 
energy E r defined by the full action S + 5S and considered as the function of the 
renormalized parameters appears to be finite both in the regularized theory, and also 
after the removing of the regularization. 

Thus, within such 'minimal addition' renormalization scheme we obtain for the 
renormalized Casimir energy the following result 

£ r = S fin + Xf r + g r X 2 r 1 (33) 

where finite parameters f r , g r must be determined with appropriate experiments, or fixed 
with the normalization conditions. The number of the required conditions is determined 
by the (in) dependence of the coupling constant A on the slab thickness £. However, in the 
most interesting cases considered below only one normalization condition is necessary. 
Indeed, let us consider the Casimir pressure. It is the observable quantity defined 

by 

V = ~ W . (34) 

For constant A it leads to 

dEfin g r X 2 

Pmatt ~ Qj- + [OO) 

§ We remind the reader that we consider a free scalar field, thus the counter-term SS is the only one 
required. 



Casimir type effects for scalar fields interacting with material slabs 



10 




where only one renormalization parameter is present. 

The original definition (j2J) of the distribution function 9{£,Xz) corresponding to 
this case can be interpreted as preserving the amount of matter within the slab: 
J dx30(£,X3) = 1. However, g r remains up to the experimental fixing since there is 
no other natural normalization condition that could be used in this case, unlike the case 
discussed in the following section. The plot of the p ma tt as a function of £ for different 
values of g r is presented at Fig. CD Note the non-monotonous dependence of the pressure 
on the width of the slab for comparatively large negative g r . 



3.5. Dirichlet limit 

Alternatively, one can consider the density of the matter to be fixed and calculate 
the pressure under this condition. Then the distribution function has a different 
normalization condition J dxs6(£, X3) = £, which is equivalent to the change of variables 
A — > £k in the formulae f l29|) -f l33|) . In this case the two counter-terms in f l33j) can be 
effectively combined into a single one g r of mass dimension one £ r = £fi n + n£g r , which 
gives for the pressure 

dSfin ~ /nr\ 

Pr = Qj- - K>9v (36) 

For fixing g r we can use what we a call a Dirichlet limit, that also provides a natural 
way to establish a correspondence between our results and the previous calculations. 

One can note that putting k = —m 2 and then taking the m — > 00 limit effectively 
converts the system under consideration into a massless scalar field confined within a 
finite interval in x^: £3 G {—1/2,1/2) subject to Dirichlet boundary conditions at the 
endpoints. 

To investigate the behavior of the finite part of the Casimir energy Sfi n in this 
limit, one has to construct carefully the m — > 00 asymptotics of ( 1291) . It can be done 
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rigourously with help of the Taylor expansion of the integrand of (12 9 p in powers of e 2mpt 
and consequent resummation of the same order contributions. This procedure yields 

m 4 £ ( 7r 4\ m 3 n 2 



^ = -lW + U"9j4^- 1440^ + ° (m " ) (37) 

Now we can require that the renormalized Casimir pressure ( 1361) in this limit coincide 
with one calculated for the case of a massless scalar field subject to the Dirichlet 
boundary conditions [9] 

= (38) 

This condition fixes the renormalization parameter g r of ( )36l) 

o 

771 

* = -128^- < 39) 
The Dirichlet limit procedure presented here looks somewhat similar to the 'large- 
mass prescription' widely used in the calculations of the Casimir energy [6]. However, 
the two schemes have different physical motivation. While in the large mass prescription 
it is argued on some general grounds that the Casimir energy must vanish in the large 
mass limit, in our case we collate a particular limit of our results with a well known 
(unambiguous) physical situation. 

3.6. Vacuum Instability 

It is well known that quantum systems might become unstable for sufficiently deep 
negative potentials. In our case the threshold is given by A = —£m 2 and for deeper 
potentials the Casimir energy ( 129]) acquires an imaginary part [26] . However the 
generation of the imaginary part is not only regulated by the strength of the potential 
but also by its width. Moreover, the stability of the systems depends on the separation 
£ in a non-monotonous way. 

From the technical point of view the imaginary part of the Casimir energy ( 1291) is 
due only to the negative values of the argument of the logarithm in S(p). It is clear 
that S(p) is still real for n/m 2 = —a > — 1. However, for bigger values of a: a > 1, 
the argument of the logarithm is not positively defined anymore in the region of small 
p: p < my/a — 1. As a function of the separation it changes sign at every i§ satisfying 

1 ( 2£V „., , 

7r — arctan . (40) 



\P\ 




E 2 - \p 



The maximum separation at which the Casimir energy is still real is given by the 
minimum of £q as a function of p. It is acquired at p = and we can conclude that the 
Casimir energy is real for all separation i less or equal to 



1 ^_ arctan 2^, (41) 



m\/a — 1 V 2 — a 

while for larger £ it is complex. The dependence of the Casimir energy on the separation 
£ for £ < £ is presented on Fig. [2J 
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Figure 2. Casimir energy £fi n (|29|) as a function of £ for m = 1, k = —am 
different as. 



2 and 



3.7. Absence of divergences 

A particulary interesting case of an unstable system is given by A = —2£m 2 (or a = 2 in 
notation of the previous subsection). It is characterized by the absence of divergencies as 
it immediately follows from f l31~|) . It is possible to check that this value of A corresponds 
to the case of vanishing second heat kernel coefficient a<± — the situation known to 
possess no divergencies [I]. Note that in this case the dependence on the auxiliary 
parameter \i also disappears. 

Using the condition elaborated in the previous section we conclude that £c as 
acquires an imaginary part for Im > it/ 2. It is pedagogical to rewrite this condition in 
CI units 



then at the r.h.s. we recognize immediately the Compton wave length (up to numerical 
factor) of the particle corresponding to the quantum field. Thus, at the separations 
(much) larger then the Compton wave length the Casimir energy becomes complex, and 
the particles are created in the potential well. Such property has a clear semi-classical 
interpretation — at the separations less then the Compton wave length, there is merely 
not enough space between the semi-infinite slabs to host any particles. 

In the opposite limit - - at short separation, - - one can derive the following 
asymptotics for the Casimir energy 



It reveals that for small separations: x < x = im ~ 0.473 the potential is attractive. At 
x the energy has an absolute maximum, the force F = —md x Sc a s vanishes (at the 0(x 2 ) 
level), and the system has an unstable equilibrium position. With larger x > x (but still 
x < 7r/2) the Casimir interaction give rise to a repulsive force on the semi-infinite slabs. 



2 mc 
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This behavior has particular similarities to one described in |4] for negative second heat 
kernel coefficient a 2 < 0. However, we remind the reader, that in our case a 2 =0. This 
analytical study is supported by numerical evaluation of the Casimir energy presented 
at the Figure [21 

Summarizing, we can say that at small distances £ < £, the gap between the slabs 
tends to shrink thus restoring the homogeneity of the space filled uniformly with matter. 
If additional energy is brought to the system be moving the slabs to the distances £ > I 
the system becomes unstable emitting particles and the gap growing. This can also be 
thought as restoring the most homogenous state of the space-time filled uniformly with 
the particles. 



4. Interaction of two finite width mirrors 

4-1. The statement of the problem 

The method developed above does not only give in a closed form the main characteristics 
of a QFT model with a single material slab, but can also be generalized for more 
complicated geometries. 

In particular, let us consider two plane slabs of thickness £\ t2 interacting over the 
distance r. In this case the action can be written as 

S = S + S def (43) 
So = ^J d A xct){x){-dl + m 2 )0(x) 

d 3 x I k 2 I dx^<j) 2 {x) + K X I dxz<t> 2 (x) J 

V J — (12— f-2 J 0,1 J 

we assume independent of £3 but corresponding generalization is possible, a± + a 2 = r. 

To obtain a formal expression for the generating functional G[J] Q with S given 
above one can proceed in one of the following ways. In a most direct approach, one 
introduces auxiliary fields ^1,2 following (jlj) separately for each of the slabs, and integrate 
consequently over 0, and then over Alternatively, one can consider the system 
with one slab as an initial 'free system', and just substitute in final formulae (jSJ) and (jSJ) 
the free field propagator D with S (1131) . With the knowledge of the explicit formulae for 
S the later approach is much easier, and the results of both of them naturally coincide 



G[J] = [DetV 2 ]~ 1/2 e^ J , (44) 

S 2L = S - k 2 (S0^j V 2 (p 2 S) (45) 

V 2 = 1 + k 2 P 2 , (46) 

P 2 = 2 S0 2 (47) 



here S = D — K^DOijQi 1 (0\D) * s the propagator of the free scalar field in the presence 
of solitary slab number 1, n iy 0j and Qi, as well as other notations used below have the 
same meaning as in Sect. [3] but related to the slab denoted by the subscript. 
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The linear (integral) operators V 2 = V 2 (:c, y), P 2 = P 2 (x,y) are defined within 
the support of the layer 2. Accordingly, in what follows we assume that the free 
spatial arguments of the Fourier transformation of the integral operators belong to 
(-a 2 - £ 2 , -a 2 )- 



4-2. Casimir energy 

The energy density per unit area of the layers is given by 

S = J ^^T^H^x^y,)]. (48) 

For the explicit calculation of the energy we follow the technique developed in Sect. 13.31 
and introduce operator J 2 

J 2 = V 2 - x - 1. (49) 

The defining property of J2 can be expressed via ( )46|) with the following operator 
equation 

J 2 + k 2 P 2 (1 + J 2 ) = 0. (50) 

According to (1471) P 2 does not depend on k 2 , which allows us to express the K 2 -derivative 
of the integrand of (j4"8~|) in the following form 

° lnV 2 = ^ = P 2 (l + J 2 ) = -^ 



OK 2 V 2 K 2 

and for the Casimir energy density we obtain in complete analogy with the case of one 
slab 

4-d r d ^ f d d - i p 



To solve ( 1501) and find the operator J 2 explicitly we first use the expression of the 
propagator S (fl3"|) to calculate P 2 (jlTj) . Assuming that slabs do not intersect r > e > 
we employ f JT4|) and derive 

P 2 (x, 2 /) = P £2 (x, 2/ ) + c 1 e^ + ^ (52) 

Cl = Mi e -2a^( 1 _ e 24 v ^ l) 

2P 

Taking this expression into account we deduce that J 2 possesses properties 
analogous of those of U ( jTOl) and satisfies the following equation 

K p J 2 = -k 2 . 

We search for J 2 as a sum of partial solution of this equation and general solution of its 
homogeneous version 

J 2 (x,y) = -K 2 V p {x,y) + Ae {x+y) ^ + B [ e {x ~ y) ^ + e {y ~ x) ^) + Ce~ {x+v) ^. 
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To deduce the coefficients A, B , C we require that (150]) is identically satisfied for all x, y 
belonging to (— a 2 — L 2 , — a 2 ), then 



A—— _^ e 2(a 2+ £ 2)v ^ (^2^ + + 



(53) 



£ = - ^4 (e 2 - £ (£ - VpJ) 2 + 2 Cl n 2 E) , 



C = - 



2^ 
K2C2 



2^ 

here £ 2 is defined through 

C" 1 = £-V a2B + 2ciK 2J E(e % v^ - 1). 
For the trace of J 2 we obviously have 



TrJ 2 



J 2 (x, x)dx 



(54) 



(55) 



a 2 — £ 2 



2£oB 



>k 2 - Ae~ 2a2 ^{l - e~ 2h ^) + Ce 2a2 ^{l - e 2h ^)) 



Using explicits for A, B, C (1531) for the energy we derive 



£ 



A-d 



-/I 



k 2 



duo 



£2 



K 2 
K2?2 



d d ~ l p 



-£2, 



2^ 

^2(2 
2p 2 



2(2^)*-! 

(1 + 2( 2 [e 2a * E (E - ^p~ 2 ) 2 + 2 Ci k 2 E}) 



(56) 



1 - e 2£2 v^) ( K2 e 2a2B + 2 Cl E(2E 2 + k 2 )) . 



This somewhat cumbersome expression can be substantially simplified since the 
integration over k 2 can be made explicitly. This also makes explicit the 1 -B- 2 symmetry 
of the layers. One obtains 

d d-ip 



&2L = P 



4-d 



■X 



X lot 



2(2n) d ~ 1 



16E 2 £x£ 2 Jpi^ 



which finally is rewritten as 
&2L — £\ + £2 ' 



d 3 p 

2{2^y- 



log 



-2Er 



J] kS(1 - e 2 ^) 



i=l,2 



(57) 



Here £\ t2 give the self-energy ( 1271) of the solitary layers 1, 2 correspondingly. The third 
term in (1ST)) represents the interaction of two layers and vanishes in the limit r — > 00. 
Taking into account the p — > 00 behavior of £j ( TT2)) one notes that the interaction term 
is UV finite, and the removal of regularization made in ( 1ST)) is indeed justified. This is in 
perfect accordance with the general considerations [36] of the finiteness of the Casimir 
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interaction between disjoint bodies. Thus, the interaction between the slabs is given by 
unambiguously finite force 

dS 2L 



— ; — 
or 

(5? 



(2vr)3 e 2& -n i=1 , 2 ^(l-e 2£ ^)- 
4-3. Interaction of a single slab with a delta-spike 

Basing on the general formulae for the Casimir interaction of two layers several limiting 
cases could be considered. 

First of all, the limit £ 2 — > with k 2 = X 2 /£ 2 , ^2 fixed, brings the interaction of a 
delta-function spike with a finite-width slab 

^-Ws&^O-t^ 1 -^)- (59) 

It is worth mentioning that in this system there are two (potentially) observable 
quantities: the self-pressure of the slab in presence of a delta-spike, and the interaction 
between two of them. According to this, we retained in (159]) the self-energy of the first 
layer S\, while discarded the self energy of delta-spike which is geometry independent 
and by no means is observable. Apart from this, we note that no trace of any 'sharp 
limit' complications [29] arise in this situation. The interaction of a finite width slab 
and an (infinitely) thin one is not influenced in any sense by the divergencies associated 
with the limit £2 — > 0. 

Now we can further consider in ( |59|) the limit l\ — > with k% = \i/£i, X± fixed, 
to prove that f|5T|) does lead to the standard result for the Casimir interaction of two 
delta-spikes separated by the distance r [6] 

f f d3 f 1 (, e~ 2Er A 1 A'2 \ 

f » = /w' os r p + p + aJ • (60) 

where again the distance independent contribution was omitted. 

The force acting between two distinct objects as function of separation is plotted 
on the Fig. [3] for the following cases: two finite-width layers, finite-width layer and 
delta-spike, two delta-spikes. The most interesting peculiarity of two slabs interaction 
is that it remains finite in the limit r — > 0, which can also be seen at the analytical level. 

4-4- Connection to the Lifshitz formula 

It is also possible to consider the opposite limit of the slabs of infinite width separated 
by the finite distance r (i.e. the interaction between two semi-infinite material slabs). 

In this case we put a\ = a 2 = r/2, £\ = £2 = £, K = £&%-, and consider £ tending to 
infinity: £ — > 00 while keeping Ki finite. Then basing on S 2 l (]57j) we derive 



d^p 



e -2Er K l K 2 



(61) 
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where the r-independent (divergent) terms are omitted. 
Taking into account that 

p li2 = E 2 + k 1>2 (62) 

we can rewrite (ED further as 

. 2Er {E - y/pl){E - 7777) 



d?p 



1-e" 



This gives for the force per unit area of the semi-infinite slabs 



(63) 



/ " dr 

r d 3 p E 

(64) 



(27r) 3 r 2/?r (E+ y /pQ(E+ y /p2) _ 1 ' 

e (b-VpT)(-b-Vp3 

This formula is somewhat similar to (3.8) in [8] which represents the contribution of 
the TE modes of the electromagnetic field to the force acting between two semi-infinite 
dielectric slabs. One notes that the correspondence becomes exact if in our model we 
introduce particular dispersion into the interaction of quantum fields with the material 
defect, namely 

4f = «£f (P) = (61,2 - (65) 



This is formally equivalent to substitution of standard dispersion relation E = \/p\ + p 1 
with the new ones, Ei = a/ t%p\ + V 1 > valid inside the corresponding slabs. The new 
parameters e can be considered as the dielectric permittivity of the matter constituting 
them. 

Since the contribution from the TM modes can be obtained by replacement [38] 

^p 2 + epl 



E=dP+p&-+ ^ ^ (66) 
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we find that the following dispersion 

<f - <f (P) = - l) f + " l) Vl (67) 

is to be introduced into the coupling constant to obtain the TM contribution to the 
dielectric slabs interaction. Introducing ( 1651) and ( 1671) into ( 1611) and summing both 
contributions we immediately recover the celebrated Lifshitz formula [H] 

^/ = -4^2" / d P° I dp 2 E(al^ 1 E + d^ x M ) (68) 



Jo 



l-TM 



(E - y/pF){E - Vp™) 



(70) 



where 



p™ = -E 2 + «™(p), p£? = £ 2 + (p). (71) 

The standard form of the Lifshitz formula in terms of Si^ and p (we use the tilde 
to distinguish the Lifshitz variables) is restored with the following substitution 




p=J^ + l, §! = Vei-l+P 2 , 52 = ^2-1+^ (72) 



4-5. The self-interaction limit 

The same limiting considered above provides an insight into the choice of 

normalization condition discussed in Sect. I3.5l and lead to a better understanding of the 
ambiguities known in the computations of Casimir self-interactions. 

One notes that £i t 2 — > 00 limit taken in ( 1431) and followed by substitution K\ = k 2 — > 
—k, m 2 — > m 2 + k effectively transforms the two-slabs system into a single-slab one. 
Indeed, after this substitution the separation between the mirrors itself converts into a 
material slab in otherwise empty space, and thus one returns to the configuration of a 
single slab of finite width described classically by (1). The same limiting procedure can 
be also performed after calculation of the quantum effects in two-slabs system. Then 
we conclude that the same physical system of a single material slab can be described 
within two different (independent) approaches: directly calculating quantum effects for 
(1), or starting first with (|43|) and performing the limiting procedure after quantum 
calculations. For the self-consistency of mathematical description of a given physical 
system one would normally require the exact coincidence of results obtained in different 
approaches. 

For the sake of comparison let us rewrite (IM]1 obtained in the infinite width limit 
substituting there K\ = k 2 — > —K, and m 2 — > m 2 + k 



T 1 _fp k 2 ^E 2 + k 

L " I (27T)3 ^2 _ e 2,vTO (jg + ^W+^f 
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On the other hand, the self-pressure calculated as r-derivative of ( 12 7p can be written as 

r Ad-\* k (K + e 2r ^ W +^(E + y/WT^) 2 ) 
= — Kg — ^ d / , 3L 1 / J , . (74) 

where we introduced explicitly the counter-term needed for renormalization of J 7 . 
By direct comparison one immediately notes that ( 173]) and ( 1741) differ by a 

A-d r 



= -K- g + ^ {m2 + K)d ' 2 2d -^ r{ - d/2) . (75) 

It is easy to check that A vanishes in the limit d — > 4 if we fix the counter-term g 
following the prescription advocated in Sect. I3.5[ and the two results obtained in two 
different ways coincide. 

Thus, we can say that two approaches describing a single-slab system differ by the 
presence of divergencies. In the direct approach based on action (1) the self-pressure 
( 1741) appears to be formally divergent, thus leading to an ambiguity represented by the 
counter-term g. At the same time, the Casimir force between the two slabs is finite as it 
generally should be for distinct bodies, [36] , and this property is preserved in the infinite 
width limit described above. This limit transforms the systems under consideration 
into a single-slab one, and thus (173]) must also be considered as a self-pressure. By 
construction it is finite and does not possess any ambiguities. For the self-consistency 
of our description we require exact coincidence of (I73|) and ( 1741 thus fixing the value 
of (otherwise arbitrary) g. Remarkably, the same value for g is obtained if we impose 
physical normalization condition based on the Dirichlet limit as discussed in Sect. 13.51 

The discussed limiting procedure establishes a new divergence free approach to 
calculate the self-pressure of a single finite width slab. However, we would like to stress 
that this limit should not be considered as a physical process but a formal mathematical 
trick. Still it allows one to obtain finite unambiguous self-pressure of a single material 
slab. 

On the other hand, one might argue that the approach based on the action (1) is 
actually describing a reacher physics as after renormalization the self-pressure contains 
an extra parameter g. Thus it could probably give answers to a wider range of physical 
questions. However, we believe that the requirement of the self-consistency should be 
used instead to fix the value of this parameter especially taking into account that it 
leads to same value of g as derived in Sect. 13.51 Moreover, the less parameters one has 
in a model the more powerful are its predictions since in this case one has less fitting 
freedom to coincide with experiment. 

The final answer discriminating between two approaches is to be given by an 
appropriate experiment. Such an experiment could be based on measuring the change 
in energy as a material slab expands or shrinks. The feasibility of such measurements 
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can possibly be considered in condensed matter physics where the self-pressure of a 
(soft) material slab can be significant for the properties of the states of the systems. 

5. Conclusion 

In this paper we constructed QFT model of the scalar field interacting with the 
bulk defects concentrated within slabs (mirrors) of finite width. Within a rigorous 
mathematical approach we constructed basic observable quantities for such systems and 
also thoroughly investigated physically interesting limiting cases. 

For the system of a single mirror we calculated explicitly the propagator, the S- 
matrix elements and the vacuum determinant (Casimir energy). The latter one coincides 
with results obtained previously in [26], [27] within a different approach. The explicit 
formulae for the full propagator and the S-matrix elements are given for the first time. 
The Casimir energy is UV divergent and for its regularization we applied dimensional 
regularization. It allowed us to extract the finite part and to construct the counter- 
terms. The renormalization procedure requires generally two normalization conditions 
to fix the values of the counter-terms with the appropriate experiments. It is shown that 
the Casimir pressure in the system can be calculated in two different ways: for fixed 
density of matter and for fixed amount of matter of the slab. In the latter case we were 
also able to establish a physically sensible normalization condition by comparing the 
system to the well known case of Dirichlet boundary conditions imposed on the quantum 
fields. We also considered the non-stable potentials, when the creation of particles out 
of vacuum can be expected, and investigated the properties of the divergence free type 
of defects. 

The method used for the single mirror system proved to be applicable to more 
complicated systems. In particular we considered the configuration of two mirrors of 
different widths separated by the interval r. Using the technique developed for the single 
mirror case we were able to calculate explicitly the Casimir energy of the quantum field 
interacting with such defect. The system possesses two different types of observables - 
the interaction of two slabs and self-interaction of one slab influenced by the presence 
of the other. The self-interaction can also be computed in two different ways mentioned 
above. Based on the general formula we considered several limiting cases. We presented 
corresponding explicit expressions and numerical evaluations of the interaction of a single 
slab with a delta spike (which further reproduces correctly the limit of two delta-spikes) 
and interaction between two semi-infinite slabs. Both results are presented for the first 
time. Investigating the final formulae we were able to establish a direct correspondence 
of obtained results with the well known Lifshitz formula which supports the validity of 
our approach. 

Moreover, investigating the Casimir interaction of two mirrors we developed a 
formal limiting procedure which recovers in an unambiguous and divergence free manner 
the Casimir self -energy of a single slab. As well known, the separation dependent part of 
the interaction between two disjoint bodies does not possess any ultraviolet divergences. 
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This property was found to be transferable to the single mirror case by the developed 
mathematical trick thus leading to well-defined finite self-energy. 

Despite it is a unique case by now, we believe that such limits should be searched 
for in other demanded cases of Casimir self-interaction, and especially for the case of 
dielectric ball. We expect that further generalization of the method proposed in this 
paper to the case of QED and its application to the spherical geometry is also possible. 
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Appendix A. Regularization dependence 

To compare different calculations in (renormalizable) quantum field theory one should 
normally operate with renormalized quantities subject to the same normalization 
condition. However, in a number of Casimir type problems (e.g. massive scalar 
field interacting with shells, etc.) the understanding of physically viable normalization 
conditions is still lacking. Yet a sound comparison is possible if the same regularization 
scheme is applied. 

When there are different regularizations used in different calculations, the 
renormalizability of considered models still guaranties the possibility of comparison. It 
is due to the fact that in renormalizable theories regularization as such is not attributed 
with any independent physical meaning and different regularization schemes must lead 
to the same final answer [jjj. Then an interchange of regularization schemes is possible and 
well justified. In a regularized quantity one should introduce additional regularization 
(which is always possible) and then remove the original one. It is a technical problem 
to the investigator to realize this programme explicitly. 

To illustrate this idea, we consider our results for the Casimir self-energy of a 
single slab. We can compare (1271) obtained in dimensional regularization with previous 
calculations of [26J and (27] where (^-function regularization was applied, and (33] where 
cut-off regularizations is used. 

In particular one can obtain the generalization of (8) [26] (written for d = 3) to 
arbitrary (spatial) dimension d in the form 

_ 1 ( r(d/2)T( S -d/2)Q d ^ ^ m 2 )d /2-s\ rAl x 
V eff--2 ds { 2T(s) 727F^>" + m) j (A ' 1} 

| We put aside an important problem of anomalies which is concerned with breaking different 
symmetries by using different regularizations. We do not discuss this issue here, as it is not relevant 
for scalar fields under consideration. 
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where Qd = 2n d ^ 2 /T(d/2) is the volume of the (d — l)-dimensional sphere in the d- 
dimensional space. Now we can remove safely the (^-regularization putting s = 0, and 
then using (15) of [26] we arrive at 

p rwf i,i '-™ ,1 1 l "'" w (A ' 2) 

For considered case of a single plane slab sn(ik) is defined in (22) [26]. Integrating by 
parts in (IA.2[) and changing integration variable k = \J q 2 + m 2 , one can check that the 
right hand sides of (l2Tj) and ( IA.2j) coincide up to change of notation L = e, Vq = r], 
d — > d — 1. In a similar way one can verify that result obtained in [27] for the plane slab 
in 1 + 1 dimensions with help of zeta-function regularization agrees with ( 12 7\i . 

Analogously we can collate our results with those presented in Section III.B of 
[33] . The authors of [33] consider (67), (68) [33] as being regularized by cutoff function 
A in fixed number of dimensions d. Equivalently it can be viewed as dimensionally 
regularized for arbitrary (non growing) A, and thus compared with our result ( J27l) . One 
can easily check that the massless limit of (12 7p does not coincide with (67) [33J. As we 
argued in [31] the Casimir energy of a single slab presented in [33] is incorrect. 

This consideration shows that when comparing results of different calculations one 
should be aware of subtleties concerned with application of different regularization 
schemes. On the other hand, it should not be considered as an insuperable obstacle 
to a comparison. 
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